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Abstract —Model-based decompositions have gained consider¬ 
able attention after the initial work of Freeman and Durden. This 
decomposition which assumes the target to be reflection symmet¬ 
ric was later relaxed in the Yamaguchi et al decomposition with 
the addition of the helix parameter. Since then many decomposi¬ 
tion have been proposed where either the scattering model was 
modifled to fit the data or the coherency matrix representing the 
second order statistics of the full polarimetric data is rotated 
to fit the scattering model. In this paper we propose to modify 
the Yamaguchi four-component decomposition (Y40) scattering 
powers using the concept of statistical information theory for 
matrices. In order to achieve this modification we propose a 
method to estimate the polarization orientation angle (OA) from 
full-polarimetric SAR images using the Bellinger distance. In 
this method, the OA is estimated by maximizing the Bellinger 
distance between the un-rotated and the rotated T 33 and the T 22 
components of the coherency matrix [T]. Then, the powers of the 
Yamaguchi four-component model-based decomposition (Y40) 
are modifled using the maximum relative stochastic distance 
between the T 33 and the T 22 components of the coherency 
matrix at the estimated OA. The results show that the overall 
double-bounce powers over rotated urban areas have significantly 
improved with the reduction of volume powers. The percentage 
of pixels with negative powers have also decreased from the 
Y40 decomposition. The proposed method is both qualitatively 
and quantitatively compared with the results obtained from the 
Y40 and the Y4R decompositions for a Radarsat-2 C-band San- 
Francisco dataset and an UAVSAR L-band Bayward dataset. 

Index Terms —Synthetic aperture radar, radar polarimetry, 
polarization orientation angle, stochastic distance 

1. Introduction 

OLARIMETRIC target decomposition is a technique 
to characterize scattering mechanisms from polarimetric 
synthetic aperture radar (SAR) data. Target decomposition 
techniques can be broadly classified into two categories: (1) co¬ 
herent decomposition techniques which utilize the information 
contained in a target scattering matrix [S], and (2) incoherent de¬ 
composition techniques which utilize the second-order statistics 
in terms of covariance matrices ([C] or [T]) derived from the 
scattering matrix [S]. Furthermore, incoherent decomposition 
techniques can be subdivided into eigenvalue/eigenvector based, 
and model-based. Eigenvalue/eigenvector based decompositions 
provide a unique solution in terms of the scattering mecha¬ 
nisms Q, 0. However, Ref. 0 raised questions regarding 
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the assignment of each eigenvector to one of the independent 
scattering mechanisms by the spectral decomposition of the 
coherency matrix. Nonetheless, the solutions provided by 
model-based decompositions depend on the assumptions made 
about the physical scattering model. 

Model-based decompositions have gained considerable at¬ 
tention after the initial work of Freeman and Durden Q . Their 
decomposition assumes the target to be reflection symmetric, 
i.e., that the co-polarized and the cross-polarized components 
are always uncorrelated. This assumption was later relaxed 
in the Yamaguchi et al decomposition (Y40) 0 which 
included the helical scattering as a fourth component. These 
two decompositions are widely used in both practice and 
literature because of their simplicity and computational ease. 
Besides these, the decomposition of Arii et al 0,0 and 
Neumann et al ® 0 introduced different scattering models 
for vegetation volume component. Jagdhuber et al 0 
proposed a full-polarimetric decomposition with multi-angular 
data acquisition to estimate soil moisture from bare ground 
and vegetated soils. Cui et al E) proposed a complete and 
a exact decomposition of the coherency matrix into a volume 
and two single scattering with non-negative scattering powers. 
Lee et al m investigated the shortcomings of model-based 
decompositions, and suggested several models to alleviate them. 
Recently, Jagdhuber et al HD developed a hybrid model-based 
and eigenvalue/eigenvector-based polarimetric decomposition 
technique with generalized volume model for soil moisture 
estimation under vegetation cover. 

A major advancement was made with the orientation 
compensation application in model-based decompositions. This 
was necessary because of the fact that a target with different 
orientations in the plane orthogonal to the radar line of sight 
(LOS) will have different polarimetric responses. A number 
of decomposition methods with orientation compensation have 
been proposed 0 0 -p8|. The fundamental idea behind such 
compensation is to minimize the cross-polarization component. 
The orientation compensation can, to a certain extent, reduce 
the overestimation of the volume power in a model-based 
decomposition and increase the double-bounce power. Most 
notably among these methods are Lee et al ||T^, the three- 
component model-based decomposition by An et al 
the orientation compensated four-component decomposition 
by Yamaguchi et al (Y4R) |T6| , and the generalized four- 
component decomposition by real and complex rotation of the 
coherency matrix by Singh et al 0 - The decomposition 
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methods of Arii et al. Q and Chen et al. HD also use 
orientation angle in their scattering models. 

In general, orientation angle (OA) estimation methods can be 
broadly categorized into two groups according to the input data: 

( 1 ) OA derived from a Digital Elevation Model (DEM), and 

( 2 ) OA derived from PolSAR data. The slope and the azimuth 
obtained from the DEM is used to estimate the OA. Apart from 
the DEM derived orientation angle there are few other methods 
available in the literature which directly use PolSAR data to 
compute the orientation angle. The phase difference between 
the RR-LL (Right-Right and Left-Left) circular polarizations 
has been used in fT^ to estimate the orientation angle. The 
polarization orientation shift is used to infer terrain slopes p0| . 
In a later study, Kimura et al. computed the shifts in the 
polarization orientation angle in built-up areas. The beta {P) 
angle obtained from the eigen-decomposition of the coherency 
matrix in Cloude-Pottier decomposition is also used as a 
measure of the orientation angle. For deterministic (coherent) 
scatterers, the orientation angle (diagonalization angle) is 
obtained from the Cameron et al. decomposition p^ . Xu et 
al. estimate the target deorientation angle by minimizing 
the cross-polarization. 

In this paper we propose to modify the Yamaguchi four- 
component decomposition (Y 40 ) scattering powers using the 
concept of statistical information theory for matrices. In order 
to achieve this modification, we first introduce the concept of 
stochastic divergence in Section |I^ In Section |I^ we propose a 
method to estimate the polarization OA from full-polarimetric 
SAR images using a stochastic distance. In this method, the 
OA is estimated by maximizing the stochastic distance between 
the un-rotated and the rotated T33 and the T22 elements of 
the coherency matrix [T]. In Section IV the powers of the 
Yamaguchi four-component model-based decomposition (Y 40 ) 
are then modified using the relative stochastic distance between 
the T33 and the T22 elements of the coherency matrix at the 
estimated OA. The results show that the overall double-bounce 
powers over rotated urban areas improve significantly with the 
reduction of volume powers. The percentage of pixels with 
negative powers is also reduced from the Y 40 decomposition. 
In Section |V| the proposed method is both qualitatively and 
quantitatively compared with the results obtained from the 
Y 40 and the Y 4 R decompositions. 


II. Stochastic Distances 

In this section we introduce our main tool: the concept 
of stochastic divergence p 4 | . Such class of techniques has 
found many application ranging from classification | [ 25 | , cluster 
analysis and goodness-of-fit tests p 7 | . Specifically in SAR 
image processing and analysis, these statistical separability 
measures have been used for estimation (281 , classification 
and segmentation 1^- m noise reduction ( 3 ^ , and change 
detection | |^ . 

Measures of divergence among probability distributions are 
an intuitive approach for comparing models. In this work, 
we use the concept of stochastic divergence to estimate the 
polarization OA from full-polarimetric SAR images. 

Goudail and Refregier ( 3 ^ used numerical integration to 
compute the Kullback and Bhattacharyya distances to derive a 


scalar measure of contrast between areas under the 2 D circular 
Gaussian model. Nascimento et al. (291 used the same approach 
for intensity data under the distribution, and compared the 
performance of test statistics derived from a number of (h-cj)) 
distances (Kullback-Leibler, Renyi, Hellinger, Bhattacharyya, 
Jensen-Shannon, Arithmetic-Geometric, Triangular, and Har¬ 
monic Mean) as features for image classification. Cintra et 
al. | [ 35 | compared some of these measures with parametric and 
nonparametric tests, and they outperformed other techniques 
in terms of efficiency and robustness. 

Frery et al. | [ 3 Q| , (361 obtained analytic expressions for en¬ 
tropies and distances between Wishart models for full PolSAR 
data. Besides being able to build adequate classificators, as the 
one presented by Silva et al. such expressions lead to 
efficient edge detectors (771 and nonlocal means filters 1^ . 

The cornerstone of these results is the family of (h-cj)) 
divergences, as defined by Salicru et al. ( 3 ^ . Consider two 
probability distributions Vi and V2 defined on the same 
support S and characterized, without loss of generality, by the 
densities fi and /2. Given any strictly increasing (decreasing, 
respectively) function h: R+ ^ [O5O0] such that h{0) = 0, 
and any convex (concave, resp.) function (j)\ R+ ^ [0, cxo] 
then 

= (1) 

is the divergence between the distributions. Mild properties 
are required for the defining functions h and 0, namely that 
lim2.^o+ exists, that 00(0/0) = 0, and that for any a > 0 
holds that limg^o+ e 0 (a/e) = alima^^oo = O 0 (a/O). 

Any divergence can be turned into a distance making 
d^(Di,D2) = (£>5(Pi,p2) + D’^{V2 ,Vi))/2 . These are 
stochastic distances as they have the properties of symmetry, 
non-negativity and identity of indiscernibles. Besides that, as 
attested by the aforementioned literature, they can be turned 
into convenient tools for image analysis. 

Adequate choices of h and 0 lead to many well known 
distances, being the Hellinger distance the one we will use 
in this work. It is obtained with h{x) = x /2 and 0(x) = 
(x^/^ — lY, then ([T]) becomes 

4(pi,p2) = i-y' vTT, 

as it is also a distance. 

In the following we will instantiate this distance under the 
most important model for ful PolSAR data: the Wishart law. 

Denote 2; = [Si^ 82^ • • - ^ SpY di complex random vector with 
p polarizations (p = 3 for reciprocal medium or monostatic 
radar) where Si, 82 and 83 are either Shh, V^Shv^ Syv 
or Shh + 8yy, Shh ~ 8yy, 28 hv the elements of 

a scattering vector in the Lexicographic and Pauli bases, 
respectively. In PolSAR data analysis, 2; is often assumed to 
obey a zero mean multivariate complex Gaussian distribution 
characterized by the following density: 

( 2 ) 

where Y) is the Hermitian positive definite covariance matrix, 
u* denotes the conjugate of v, and its transpose. In order to 
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increase the signal-to-noise ratio, L independent and identically 
distributed samples are averaged to form the L-looks covariance 
matrix 


definite coherency matrix [T] is unitarily rotated by [U3]: 
TO] = [U 3 ][T][U 3 ]-^ 
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( 3 ) 
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The p X p Hermitian positive definite matrix Z follows a 
scaled complex Wishart distribution whose density is 

^ I>--|2|^--’exp[-LTr(S-Z)] 

r,(L)|s|‘ 

where Tr is the trace and Tp{L) = Ylu=i +1). 

The Hellinger distance between two scaled complex Wishart 
distributions with the same number of looks Li = Lj = L is 
given by 


[U 3 ] = 


( 7 ) 


1 0 0 
0 cos( 26 >) sin( 26 >) 

0 — sin( 26 >) cos( 26 >) 


The orientation angle 0 G [— 7 r/ 8 , 7 r/ 8 ] is estimated by 
minimizing the Tss{ 0 ) component, i.e., dT^^{ 0 )/d 0 = 0, with 
this: 


2 733 — T22 


( 8 ) 


The effect of orientation on the three diagonal elements of the 
coherency matrix is as follows: 




1 - 


2\/|^i||^il 


L 


( 5 ) 


where we indexed the distributions by their parameters. Con¬ 
sequently the Hellinger distance for the 1 -dimensional (p = 1 ) 
multilook intensity distribution is given by 


dH{crh<^l) = 1 - 


( 7 i( 7 i 


2K2 + <t2) 


( 6 ) 


Expression is the Hellinger distance between two Gamma 
distributions with same shape parameter L and expected values 
erf and cr|. In the context of PolSAR data we can associate 
cr? and of the two probability density functions with the 
diagonal elements Tu, Tu{ 0 ) for ^ = 1, 2 ,3 of the un-rotated 
and the rotated Hermitian positive definite coherency matrices 
respectively; cf. 0. 


HI. Polarization Orientation Angle Estimation 


1 ) Til = |HH + VV |^/2 is roll invariant for any 0 , 

2 ) T22 = |HH — VV |^/2 always increases or remains the 
same after the OA compensation, and 

3 ) T33 = 2 |HV|^ always decreases or remains the same 
after OA compensation. 

In this work, we have used a statistical information theoretic 
measure to estimate the OA from full-polarimetric SAR 
imagery: the OA is estimated by maximizing the Hellinger 
distance between the T33 and the T22 elements of [T]. 

The proposed method estimates the OA by first maximizing 
the Hellinger distance between the un-rotated (al, erf) and the 
rotated (crK^), elements, over the [— 7 r/ 4 , 7 r/ 4 ] range, 

leading to two candidate angles: 


03 = argmax ^ 1 — 

— 7r/4<0<7r/4 


02 = argmax ^ 1 — 

— 7r/4<0<7r/4 


2ycr|cr|(6>) 

3 

_al + aUe)_ 

J 


3 

erf + (7|(6») 

J 


and 


( 9 ) 


In this section we estimate the polarization OA from full- 
polarimetric SAR data by using a stochastic distance measure 
between the elements of the coherency matrix, which is 
assumed to follow a complex Wishart distribution. 

The OA is zero for reflection symmetric media, but OA 
shifts are induced for surfaces with azimuthal tilts and for 
buildings oriented perpendicular to the radar LOS. Apart from 
these, the OA is appreciable for low frequency (L and P 
band) data in forested areas due to surface topography. In 
this context, a recent adaptive-model based decomposition 
with topographic polarization orientation compensation (TPOC) 
has been proposed in p 9 | . In this, the volume orientation is 
removed from the generalized volume component, and the 
conventional adaptive model-based decomposition is modified 
by introducing the TPOC concept. 

In general, the primary effect of compensating the OA is in 
the reduction of the cross-polarization (T33) component and 
the increase in the co-polarized (T22) component. Eollowing 
the proposal of Lee and Ainsworth 0. the Hermitian positive 


Two maxima are found at 0 = 0{3,2} and 0 = 0{3,2} i 7r/4, 
and the OA is chosen such that the Hellinger distance (dnepa) 
corresponding to is greater than the Hellinger distance 
(dHcf)2) corresponding to either at 0 = 0{3,2} or at 0 = 
0|3 2} ± 7r/4. This condition corresponds exactly to the case 
mentioned earlier, where the cross-polarized component is 
minimized, whereas the other peak corresponds to the situation 
where the cross-polarized component is incorrectly maximized. 
Einally, the OA (Oq) is obtained by wrapping 0 in the range 
[— 7 r/ 8 , 7 r/ 8 ] using which is then compared with the OA 
obtained by the method stated in 0: 


Oo 


0 + 7r/4 if 0 < —7r/8, 
<0 — 7r/4 if0>7r/8, 

0 Otherwise. 


( 10 ) 


The implementation of the proposed method is given in 
Algorithmic The code is freely available at https://github.com/ 
avikesre/SD Y 40 / 
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(a) (b) (c) 


Fig. 1. (a). UAVSAR L-band Pauli RGB image of the Hayward area, (b). OA estimated by the method proposed in [l4\ and (c). OA estimated by the proposed 
method using the Hellinger distance. 


Algorithm 1: Orientation Angle Estimation 

Data: [T] matrix 

Result: Orientation compensating angle (6>o) 

1 begin 

2 


r]' ^ [U3(0)][T][U3(0)]-' [-7r/4,+7r/4] ; 

3 

Find common peaks in ^ dH{cF‘i^cF‘i{0)) and 


d 

H(f)2 ^ dH{o-2,CF2{d)) ; 

4 

At 0{3,2}- if dHcf)3 > then 

5 


Select (p ^ (ps ; 

6 

else 

7 

1^ Select (j) ^ (j) 2 \ 

8 

if 0 > 7r/8 then 

9 


6>o ^ 0 - 7r/4 ; 

10 

else if 0 < —7r/8 then 

11 


6*0 •«-?!’ + 7 r /4 ; 

12 

else 

13 

|_ do ^ (p ; 





Fig. 1^ to estimate the OA. It can be seen that the OA is also 
correctly estimated by using the KL distance for the coherency 
matrix given in ( pT] ). Our choice was based on computational 
cost, since the Hellinger distance is the least expensive among 
the available ones. 



Fig. 2. Variation of Hellinger distances (d/f (fj?, (To (0)) in red and 
blue) and KL distances {dKL{(^,(^{4>)) in cyan and 
dKL{(^ 2 T^ 2 ^^)) i^ green) with (j)\ 6 h is the relative Hellinger distance 
between T 33 and T22 at cf) = 60. 


In the following we use a coherency matrix [T] extracted 
from an urban area, to illustrate the proposed methodology: 


[T] 


4.56 2.28 + 0.72i 0.02 + 0.67i 

2.28 - 0.72i 6.06 1.90 + 0.27i 

0.02 -0.67i 1.90 -0.27i 3.50 


( 11 ) 

The maxima estimated by the proposed methodology are, 
approximately, at (j) = —31° and at 0 = 14°, as shown in 
Fig. The orientation compensating angle at Oq = 14° is 
chosen because it satisfies the aforementioned criterion of 
> d/f (erf, cr| (?!■)). 

Apart from the Hellinger distance, any other distances 
mentioned in Section |n| can be used to estimate the OA. An 
example is shown with the Kullback-Leibler (KL) distance in 


IV. SD-Y40 

Presently, all the decomposition techniques that account 
for the overestimation of the volume power due to target 
orientation -|[T^ do so either by rotating the coherency 

matrix or by considering different volume scattering models. 

In this work, a modification of the scattered powers for the 
Yamaguchi four component decomposition (Y40) is obtained 
using the relative Hellinger distance. Fig. shows the relative 
distance 6h between the two Hellinger distances dni^cF^^ 
and dH{cF 2 ^^ 2 {^)) computed at (j) (the unwrapped orientation 
angle in the [—7r/4,7r/4] range). This relative distance is used 
to modify the Y40 scattering powers. 

According to Lee and Ainsworth |T4] |, the amount of increase 
in double-bounce power is not equal to the amount of decrease 
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in the volume power. It has been shown that the amount of 
decrease in the volume power component is usually greater than 
the amount of increase in the double-bounce power component. 
Moreover, the increase in surface power is smaller compared to 
the increase in double-bounce power. Due to the roll-invariant 
property, the helix power remains unaltered after rotation. 

In the following we quantify the effect of the proposed 
rotation on the decomposition powers of a particular target: 

0<Sh= dHial,al{(j))) - dH{(jl,crl{(f))) < 1. (12) 

From the previous section, it can be observed that the OA 
estimation is independent of L, but the relative distance 5h 
increases with increasing L until a maximum is attained, and 
then it rapidly decreases to zero as shown in Fig. 13 The Y40 
decomposition powers are then modified with the estimated 
^ where is the maximum attainable 

divergence between as and (72. 



a and P parameters for the UAVSAR L-band image are shown 
in Fig. 13 



(a) a (b) ^ 


Fig. 4. The a and the /3 parameters maps for the UAVSAR L-band Hayward 
image. 


It should be noticed that the linear mapping of the estimated 
OA to the a parameter does not guarantee maximum (and 
Ps) powers to be observed at/near 45°. The Pd and Ps powers 
are modified by using the relative Hellinger distance 6 ^, which 
is a non-linear function of the OA angle, multiplied by the a 
parameter. Figs. 5b and 5c show the 2d scatter plots of a6'^, 
and /36^ versus 0 , respectively, with the contours representing 
the density for oriented urban areas as shown in Fig. ^ It 
can be seen that the maximum density of points lies between 
(9 = 20° to (9 = 35°. 


Fig. 3. The variation of Sh with L for cj) = 14°. A maximum of 6h = 
is obtained for a certain value of L = Lm- 

In the proposed methodology (SD-Y40), we modify the 
Yamaguchi four component decomposition powers based on 
the maximum relative Hellinger distance by: 


P,"=P,(l-(5^), (13a) 

P^ = Pd PaPJ^, (13b) 

P^ = PsPPPJ^, and (13c) 

= Pc. (13d) 


where P^, Pd, Ps and Pc are the volume, double-bounce, 
surface and helix power components from the Yamaguchi four- 
component decomposition (Y40), respectively. The new four- 
component powers: P^, PJ^ and P^ are obtained by either 
deducting or adding the volume power component (P^) of 
the Y40 decomposition, adjusted by the relative distance 
respectively. 

The modulating factor 5^ is adjusted by the positive 
parameters a and P with a-\-P = 1 for PJ^ and P^ respectively 
in order to conserve the total power (TP = P^ + P^ + P^). 
Moreover, adhering to the criteria in fT4] |, i.e., the increase in 
the double-bounce scattering power is more than the increase in 
the surface scattering power, the ranges of the two parameters 
are set to 0.5 < a < 1 and 0 < /3 < 0.5. The a parameter 
for each pixel is obtained by linearly mapping the estimated 
orientation angle (0 G [0,7r/4]) to [0.5,1.0]. The maps of the 



(b) (c) 


Fig. 5. (a) The Pauli RGB of a oriented urban area, (b) and (c) the scatterer 
plot along with density showing the variation of (x5^ and with 6 

In the following section the SD-Y40 method is quantitatively 
compared with the Y40 and Y4R decomposition scattering 
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Fig. 6. Pauli RGB of (a) Radarsat-2 data over San Francisco and (b) UAVSAR 
data over Hayward. Zoomed optical images are shown alongside for some 
areas. Courtesy: Google earth. 


Y40 to 9% for SD-Y40, as shown in Fig. |^d)-(f). However, 
as expected, the powers over area “B” remain unchanged. The 
double-bounce power over the area “C” has increased from 
74% for Y40 to 82% for SD-Y40. A small change of 3% for 
the surface power from Y40 to SD-Y40 can also be observed 
over this area. This could be attributed to the water body partly 
present in the area marked as “C”. 

Also three samples from the L-band UAVSAR image of the 
Hayward region of San Francisco are analyzed. In this region, 
there are certain urban areas which are approximately 10°-20° 
rotated away from the LOS, as can be seen in Fig. [2a)-(c). 
A small sample, identified as “A”, is used from this oriented 
urban area to compare the powers obtained from Y40, Y4R and 
SD-Y40. A major improvement can be seen in “A” with the 
reduction of the volume power from 60% for Y40 to 27% for 
SD-Y40. The double-bounce power has increased from 30% 
to 52% compared to Y4R, which has only increased to 36%, 
as shown in Fig. [^d)-(f) and in Table These improvements 
can be correctly attributed to the fact that the area under 
consideration is a dense urban area rotated about the LOS. 

The surface scattering power has also improved from 9.1% 
for Y40 to 19.8% for SD-Y40. However, it can be observed 
that for the urban area “B’,’ which is facing towards the LOS, 
the scattering powers obtained from the three methods (Y40, 
Y4R and SD-Y40) are comparable, as shown in Fig [^g)- 
(i) and in Table |I^ This suggests that the proposed method 
is useful for extracting proper scattering powers in oriented 
urban areas. In the case of non-oriented urban areas, the three 
methodologies provide similar results. Another patch, marked 
as “C” in this image, is a forested area. The powers obtained by 
the three decompositions are similar with the volume scattering 
power being the dominant type. 


powers for a Radarsat-2 C-band and a UAVSAR L-band. 

V. Results and Discussion 

A fine-beam quad polarization (FQ9) Radarsat-2 C-band 
image over San Francisco with a spatial resolution of 8 m x 8 m, 
and a UAVSAR L-band data with a spatial resolution of 
12m X 12 m (multilooked 3 times in both range and azimuth) 
are used in this study. The P^, Pd and Py powers are obtained 
from the Y40 and Y4R decomposition algorithms given 
in 0 and GD, respectively. The P^, Pd and the Py power 
components corresponding to the two decompositions (Y40 
and Y4R) and the proposed modification (SD-Y40) are shown 
in Fig. for the Radarsat-2 C-band San-Francisco image. The 
Pauli RGB compositions of the Radarsat-2 and the UAVSAR 
images are shown in Fig. along with high resolution optical 
images of patches of urban areas (both aligned and rotated 
with respect to the LOS), and of a forested area. 

Three types of land cover classes are considered in the 
Radarsat-2 image: (1) Area “A” is a highly oriented urban 
region, (2) Area “B” is a region with vegetation cover, and 
(3) Area “C” is the Golden Gate bridge. Over area “A”, the 
double-bounce power has increased from 2% for Y40 to 17% 
for SD-Y40, and the surface power has increased from 3% for 


TABLE I 

Average decomposition powers over a rotated urban area (Area 

’A’ MARKED IN FIG.[^A)-(C)) 


Methods 

Ps 

Pd 

Pv 

Y40 

0.05 

0.17 

0.33 

Y4R 

0.09 

0.20 

0.26 

SD-Y40 

O.II 

0.30 

0.14 


TABLE II 

Average decomposition powers over an un-rotated urban area 
(Area ’B’ marked in Fig.[^a)-(c)) 


Methods 

Ps 

Pd 

Pv 

Y40 

0.37 

1.67 

0.16 

Y4R 

0.35 

1.69 

0.14 

SD-Y40 

0.38 

1.70 

O.IO 


A major concern associated with model-based decomposi¬ 
tions is the occurrence of negative powers in double-bounce 
and surface scatterings. Many ad-hoc techniques have been 
proposed which force the negative powers to be positive. The 
negative powers issue was discussed in | [T5| , and in | [4Q| the 
constraint of nonnegative eigenvalue was proposed to mitigate 
the problem. However, in this work, the double-bounce and 









7 



(a) Y40 


(b) Y4R 


(c) SD-Y40 



(d) Y40 Area A 



(e) Y4R Area A 


M 

Ps 

m 


m 


Pv-7^ 

L0% j 



r 


(f) SD-Y40 Area A 



(g) Y40 Area B 



(h) Y4R Area B 



(i) SD-Y40 Area B 




(j) Y40 Area C 


(k) Y4R Area C 


(1) SD-Y40 Area C 


Fig. 7. Radarsat-2 C-band San-Francisco image (a)-(d)-(g)-(j): Yamaguchi four component decomposition without rotation of the coherency matrix (Y40) for 
area ’A, ’B’ and ’C’ respectively, (b)-(e)-(h)-(k): Yamaguchi four component decomposition with rotation of the coherency matrix (Y4R) for area ’A, ’B’ and 
’C’ respectively, (c)-(f)-(i)-(l): Y40 decomposition scattering powers modified with the Hellinger distance (SD-Y40) for area ’A, ’B’ and ’C’ respectively. 


the surface powers obtained from the original Yamaguchi four 
component decomposition which show negative powers are 
modified by the two positive quantities, aPy6^ and 
Hence, it can be expected that the number of pixels with small 
negative powers in double-bounce and surface scattering can be 
made positive with the addition of the two positive quantities, 
but the number of pixels with large negative powers may still 


remain negative. A comparison of the negative powers for the 
two datasets is shown in Table [Ivl It can be seen that the 
percentage of pixels with negative powers have decreased from 
Y40 to SD-Y40 for both the datasets and is comparable to 
Y4R. 




(a) Y40 


(b) Y4R 


(c) SD-Y40 





(d) Y40 Area A 


(e) Y4R Area A 


(f) SD-Y40 Area A 




¥Ps :17.0% 

/Py :7.0% ' 




(g) Y40 Area B 



(h) Y4R Area B 



(i) SD-Y40 Area B 



(j) Y40 Area C 


(k) Y4R Area C 






(1) SD-Y40 Area C 


Fig. 8. UAVSAR L-band Hayward area image (a)-(d)-(g)-(j): Yamaguchi four component decomposition without rotation of the coherency matrix (Y40) for 
area ’A, ’B’ and ’C’ respectively, (b)-(e)-(h)-(k): Yamaguchi four component decomposition with rotation of the coherency matrix (Y4R) for area ’A, ’B’ and 
’C’ respectively, (c)-(f)-(i)-(l): Y40 decomposition scattering powers modified with the Hellinger distance (SD-Y40) for area ’A’, ’B’ and ’C’ respectively. 
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TABLE III 

Average decomposition powers over a eorested area (Area “C” 

IN FlG.[g:A)-(C)) 


Methods 

Ps 

Pd 

Pv 

Y40 

0.01 

0.02 

0.38 

Y4R 

0.02 

0.02 

0.37 

SD-Y40 

0.02 

0.02 

0.37 


TABLE IV 

Percentage oe pixels with negative powers 


Dataset 

Method 

% of Pixels with 
Negative Powers 

UAVSAR 

Y40 

17 


Y4R 

14 


SD-Y40 

13 

Radarsat-2 

Y40 

8 


Y4R 

6 


SD-Y40 

6 


VI. Conclusions 

The over-estimation of the volume power and, consequently, 
the underestimation of the double bounce and the surface 
powers in the Y40 decomposition in rotated urban areas is 
of major concern. In order to alleviate this issue, the Y4R 
decomposition was proposed. Several other decomposition 
techniques have been recently proposed which address this 
issue by using different scattering models. In this work we 
propose a stochastic distance based measure to modify the 
powers estimated from the Y40 decomposition for oriented 
targets. 

We estimate the orientation angle of the target using a criteria 
based on the Hellinger distance between the T 33 and T 22 
elements of the coherence matrix [T]. Using this stochastic 
distance, we have proposed a method that systematically 
modifies the Pg, Pd and the Py powers to obtain better 
estimates. Thus, the physical property of the target, i.e., the 
orientation angle, is used to make modifications to the powers. 
The results obtained by the proposed method are encouraging. 

An L-band UAVSAR dataset over Hayward is used in this 
work due to the presence of rotated urban areas. The dataset is 
decomposed with the Y4R, Y40 and the proposed method. It 
can be seen that, on oriented urban areas, there is an increase 
in the Pg power by 11 % and an increase in the Pd power by 
22% from Y40 to SD-Y40 with a corresponding reduction of 
the volume power. A similar analysis is also performed on a 
Radatsat 2 C-band San Francisco dataset. An increase of 15% 
in the Pd power is observed for the proposed method. However, 
for the area with forest cover, the three powers from all the 
decompositions are almost identical. The authors believe that 
due to its simplicity, the proposed method can be easily used 
to modify the Y40 decomposition output, which is among the 
most popular model-based decomposition techniques in SAR 
polarimetry. 
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